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Abstract 



neutrinos. 



fN| ■ Adding singlet neutrinos to the standard model spectrum in general gives rise to 

Z-induced flavor-changing neutral currents. We study the impact of these currents 
on matter-induced neutrino oscillations in the sun and in supernovae. While the 

C*~) ' effects for solar neutrinos are negligible, dramatic effects are possible for supernova 

o 

oo 

ON 

^ ! 1 Introduction 

a: 

Neutrino oscillations []IJ can provide an appealing solution to the solar neutrino problem, 
which seems to be increasingly difficult to explain otherwise. All four solar neutrino exper- 
iments [§, H, U, D observe a deficit of z/ e 's compared to the Standard Solar Model (SSM) || 
expectation. In particular, the apparently strongest depletion of intermediate energy neu- 
trinos makes it very difficult if not impossible to account for these observations in any 
reasonable solar model [/j. Moreover the recent helioseismological data confirms the SSM 
predictions for the solar density profile || strengthening the belief that the solar neutrino 
problem is due to non-standard neutrino properties. 

The features that are required to allow for neutrino oscillations are a non-vanishing 
mass-squared difference Ay = m? — m| (implying that at least one neutrino state is 
massive) and mixing (i.e. the neutrino interaction eigenstates do not coincide with the 
mass eigenstates). The solar neutrino problem can then be solved by vacuum or matter- 
enhanced neutrino oscillations |J with 

AZi - 10~ 10 eV 2 or A™f ~ f(T 5 eV 2 . (f) 



The upper bounds [10] on the light neutrino masses, 
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m Ve < 15 eV, m Ufl < 0.17 MeV, m„ T < 24 MeV, 



(2) 



are not in conflict with such a solution, although one has to explain why the postulated 
neutrino masses are so much smaller than the charged fermion masses. 



The most popular solution is the Majorana see-saw mechanism [II |. It requires addi- 



tional right-handed SU(2)l singlet fermions with large Majorana masses which mix with 
the Standard Model (SM) neutrinos via SU(2)l breaking Dirac masses. In its most com- 
mon version this leads to a 6 x 6 mass matrix of the form 

The entries of the Dirac mass matrix M d are typically comparable to the masses of 
the charged leptons while the entries of the Majorana mass matrix Mr do not break 
weak SU(2)l and are therefore related to some very high or intermediate New Physics 
(NP) scale A. This leads to three ultra-light neutrinos with masses of order m 2 D /mn 
(mx = (detMx) 1 ^ 3 ), and singlet admixtures in the light mass eigenstates of order m^/m^. 
There are variants of the above see-saw mechanism in which expanded neutrino singlet 
sectors lead to much larger mixing with the known neutrinos. We return to a discussion 
of these scenarios later in the paper. 

The point we wish to make is the following: In order to explain the solar neutrino 



problem (and other experimental anomalies like the atmospheric neutrino problem [12 
and the LSND results |L3], [14|) one requires small, but non-zero neutrino masses. The 
most attractive scenarios for obtaining such tiny masses involve see-saw suppression. The 
important ingredient is the presence of heavy SU{2)i J singlet neutrinos, which is rather 
common in NP models. However, mixing of the known neutrinos with singlets may also 
influence the neutrino oscillations that were invoked to solve the experimental anomalies 
in the first place, since they modify the neutral current interactions in a way that is flavor 
dependent. While the existence of this effect is generic in models that allow for neutrino 



masses, it has not been extensively investigated in the literature fllS] . It is the purpose of 
this work to give a quantitative analysis of its importance for solar and supernova neutrinos. 
We find that phenomenological constraints on singlet mixing make its significance for solar 
neutrinos marginal at best, but allow for dramatic effects in the case of supernova neutrinos. 

The paper is structured as follows: In Section |2] we introduce the formalism of Z- induced 
Flavor Changing Neutral Currents (FCNCs) and work out the changes to matter-enhanced 
neutrino oscillations in the presence of non-sequential neutrinos. Phenomenological con- 
straints on neutrino-singlet mixing are discussed in Section [| The implications for the 
Mikheyev-Smirnov-Wolfenstein (MSW) solution of the solar neutrino problem are stud- 
ied in Section f| and for oscillations of supernova neutrinos in Section || We conclude in 
Section ||. 
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2 Formalism 



2.1 Z-mediated FCNCs 

Consider a model where the lepton sector of the SM is extended in a non-sequential way 
(we first consider additional SU(2)l singlets, then briefly discuss the case of additional 
triplets in Section |2T3| ). We group the known (K) and the new (N) interaction eigenstates 
in the vector SE^ = (*&k, *&n) t , which is related to the corresponding vector of light (L) 
and heavy (H) mass eigenstates \l/ = (SJ'l, ^ h) T by a unitary transformation 

(*k) =U (*l) where U = ? T KB ) . (4) 

\v N J \^ H J \u nl u nh J v ' 

The submatrices U kl and U nl describe the overlap of the light eigenstates with the known 
interaction states and the new interaction states, respectively. 

It has been known for a long time |16| that the presence of new interaction eigenstates 
in general leads to FCNC neutrino interactions. To make this statement more precise 
consider the Neutral Current (NC) operator 

N c=\*k)(*k\, (5) 

where the sum over K — e,fi, r is implicit. Effective four-Fermi couplings for Z-mediated 
neutrino scattering are proportional to O^c- Trivially, the matrix elements of this operator 
in the basis of known interaction eigenstates are 

{^ k \O nc \^>k') = Skk', (6) 
while in the basis of light mass eigenstates they are 

(#i|0*c|*j') =(*i|*,c><*jc|tfj') = U Kl Uy. (7) 

The latter are in general not equal to 5u>, since the submatrix U kl is not unitary. Thus 
there are "flavor" -changing NCs between the light mass eigenstates. 

In most cases of interest the heavy mass eigenstates are not kinematically accessible. 
Hence we need to project the interaction eigenstates onto the subspace of light "propagat- 
ing" mass eigenstates, if we want to calculate the effective operators at low energies. We 
denote the projected states by 

\VZ) = \%)Wk) = U ki \9 1 ). (8) 
Note that these states are not orthonormal to each other, since 

<#£|tt£> =<ttd*i)<*j|*jc) = uyj m . (9) 
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In the basis of these "light" interaction eigenstates the operator O^c is represented by the 
matrix 



(10) 



which is not in general diagonal since the basis {^1} does not span the space on which 
Onc is defined. 

As a simple example consider the case of two known interaction eigenstates v e and 
(I = pi or r) and one new SU(2) L singlet u s . The vector St 7 = (v e , is e , v s ) T is connected 
to the vector of mass eigenstates * M = u 2 , Vh) T ■, where v x and v 2 are light while v h is 
heavy, by a unitary transformation 









-[ 







u el u e2 u eh \ 

U a U a U eh 
Usi U S2 U Sh 




(11) 



The "propagating" neutrinos that are produced in low-energy charged- current interactions 
together with the charged leptons e and £ are 



Uel U e2 



(12) 



i.e. we have projected v e and vi onto the v\ — v 2 plane. Note that since U kl is not unitary, 
and vf are not orthogonal to each other, 



(<wn = u* el u n + u* e2 u e2 = -u* eh u eh: 



(13) 



and are not properly normalized. Thus we cannot use these states as a proper basis for the 
description of neutrino oscillations. Instead, we choose an orthonormal basis \ vf)} 

where 



_ U el \vi) + U e2 \v 2 ) 

re I 



el\ 2 +\Ue2? 



(14) 



is aligned with \vf) and 



-u^±uM 

^\Uel\ 2 +\U e2 \> 



(15) 



is orthogonal to \ In this basis the off-diagonal matrix elements of O^c have magnitude 



v°\OncW) 



0\ 



\(U* el U ei + U: 2 U e2 )(-U n U* e2 + U&U t 



el J 



\u el \ 2 + \u e2 \ 2 



\uiu 



eh u m 



(16) 
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The approximation uses (|T3|) and the fact that experimental constraints, that will be 
discussed in Section |3], imply \U e h\, |^si|> \Usz\ ^ I^SfiJ — 1? from which it follows 
that 

/t/el ^ e2 



det£7 



det . 



l^hl, (17) 



or | — UnU e2 + U£2U e i\ ~ 1. Note that in the oscillation-basis {^} the operator O^c also 
has non-universal diagonal couplings 



ln°\n \iP\ - 1 1 77 1 2 | \£ehfyh\ 
\V e \V NC \v e ) — L - \U eh \ + - — — 

1 — \Ueh 

{V t \0 N c\v t ) = : 



2 

12 ; 



U. 



eh 
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The off-diagonal element in ([T6j) and the deviations from unity of the diagonal matrix 
elements in flT%|) are small since they are quadratic in the mixing between the doublet and 
singlet neutrinos. 

2.2 NP-coupling of the W 

The couplings of the H^-boson are also affected by the presence of a heavy singlet neutrino. 
The effective four-Fermi couplings for VT-mediated neutrino scattering in the sun or a 
supernova are proportional to the Charged Current (CC) operator 

C C=We)(v e \. (19) 

In our prefered basis {v®,^®} for the description of neutrino oscillations its only non- 
vanishing matrix element is 

("?\0CCW?) =(v?We)(VeW?) = 1 - \U eh \ 2 . (20) 

The other matrix elements vanish since vf is orthogonal to v e by definition. Thus also in 
our orthonormal basis Oqc only projects onto the electron neutrino, but with a prefactor 
that is slightly smaller than the standard one. 

2.3 Adding Triplets 

Additional SU(2)l singlets are not the only relevant extension of the lepton sector. We 
briefly investigate the effects of neutrino mixing with the neutral component of an SU{2) l 
triplet. Fermionic triplets arise naturally in the context of supersymmetric extensions of 
the SM. For example, the SU(2)l gauginos form such a triplet with hypercharge Y — 0. 
However, the neutral component has weak isospin t%L = so it does not couple to the Z 
and the resulting formalism is the same as for singlets. 
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The neutral component of a triplet with Y = ±1 does couple to the Z. A good 
example is the superpartner of the neutral component of the Higgs triplet (with Y — — 1) 
in supersymmetric left-right symmetric models. It has t^i = 1 so we have to replace the 
NC-operator defined in ([5]) by 

0' NC = NC + 2\u N )(u N \. (21) 

In the context of our simple example with two ordinary neutrinos and one exotic, the 
off-diagonal element of 0' NC in the oscillation-basis is given by 



CW \,,0\ V,, l,,0\ (Uel^Nl + N2){-U*NiUt2 + ^7V2^el) / 00 \ 
U NC\ U t / -\ v e Wn)\^nWi ) - lTT ,n : lTT | 9 , {^^j 



\u el \ 2 + \u e2l - 

which is quadratic in the mixing between ordinary and exotic neutrinos, as in 



2.4 Effective Z-induced NP-couplings 

In this section we compute the effective four- Fermi couplings, G N , for the neutrino flavor- 
changing reactions uff — > vf f (/ = e, u, d) mediated by the Z, in the presence of a 
neutrino singlet. We note that all three couplings are determined by the single parameter 

e={v°\0 NC \v°), (23) 



which essentially gives the ratio between the Z-induced flavor-changing NC amplitudes 
and the usual flavor diagonal NC amplitudes. For simplicity, we will ignore the possibility 
of CP violation, taking e to be real throughout. The flavor- changing four- Fermi couplings 
follow directly from known results |]T7| for the flavor- diagonal NCs, which are reviewed 
below. 

The potential Vnc due to Z-mediated neutrino scattering off the thermal background 
fermions can be deduced from the relevant (Fierz rearranged) four-Fermi interaction 

Hf = ~ lB9 f A )f P7*(l - 75>- (24) 

The vector and axial couplings are 



g f v = 4 - 2qf sin 2 6 W (25) 
9a = 4l, (26) 

where t^ L is the weak isospin of the fermion / and qf is its electro-magnetic charge. We as- 
sume unpolarized background fermions and therefore only the 70 component of the fermion 



6 



density can contribute. The 7075 mixes "small" and "large" components of the spinor, so 
the axial coupling does not contribute. For the large component 70 is 1, thus the charged 
fermion coupling reduces to 

IlM - 159a) f = <WV/(4 - 2qf sin 2 6 W ), (27) 

where Nf is the fermion number density. Hence the term describing NC interactions in the 
effective Hamiltonian for neutrino propagation is 

C N 

H f NC = -^(2& - V sin 2 6 W ) x \v K ){v K \. (28) 

Note that the v r f , - = ^-v coupling yields the unit matrix \vk)^ v k\ hi the basis of weak 
neutrino eigenstates, which is just Onc- The factor V^ c multiplying Onc takes the values 



V NC 


= G F N e (4:sm 2 e w - 1)/V2 




V NC 


= G F N u (~8/3sm 2 e w + l) i 


/A 


v d 

V NC 


= G F N d (A/3sm 2 6 w - l)/\ 


/2, 




= G F A^(-4sin 2 ^ + l)/v 


% 


V NC 


= -G F N n /V2. 





(29) 
(30) 
(31) 
(32) 
(33) 

The contributions for the nucleons are obtained by summing the quark potentials, i.e. 
V* c = 2V% C + V$ c and V$ c = V% c + 2V# . Note that the electron and the proton 
contributions cancel in electrically neutral media. 

In the presence of neutrino singlets, the effective flavor-changing four-Fermi couplings 
G^ are given by 

rf _ (y?\n f NC \v?) _ vie 

° N = V2N f ~ VW f £ - (34) 
Equivalently, the flavor- changing NC potentials are 

Vfcnc = £ Vnc- (35) 

Note that the effective flavor-diagonal NC potentials are decreased slightly by the deviations 
of the flavor-diagonal matrix elements (v°\OncW°) an d (^I^ncW?) in eq. (p~8[ ) from 
unity. 
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2.5 Modifications to Matter Oscillations 



The MSW mechanism |§ provides an elegant solution to the solar neutrino problem and 
it might also be important for neutrinos that are produced in a supernova explosion. In 



Ref. |L8[ modifications of the MSW solution due to flavor- changing neutrino interactions 



were considered. Although the new non-diagonal four-Fermi couplings were suggested to 



stem from the exchange of new heavy particles, the results obtained in |18| can also be used 



to discuss matter-enhanced neutrino oscillations in the presence of Z-induced FCNCs. 

We have shown that the projection of the known interaction eigenstates onto the sub- 
space spanned by the light mass eigenstates results in flavor-changing NC neutrino inter- 
actions of strength e with respect to the flavor-diagonal NCs. Thus, knowing the standard 
NC-contributions it is straightforward to obtain the correct expressions for the off-diagonal 
terms in the effective Hamiltonian for neutrino propagation in matter 

z/°\ 1 /-Acos2^ + A' ASH120 + 5 \ ( v° 



Un \v°)- AE { Asin20 + B Acos20-A'J \Jf)' (36) 

Here E is the neutrino energy, A = m\ — m\ is the mass-squared difference of the two 
vacuum mass eigenstates, 9 is the vacuum mixing angle, and 



A' = A 1(1 



\Ueh?) 



2N P 



\Ueh\ 



\u:„u. 



h u ih\ 



i - u, 



eh 



Ue2Ull + UelUl2\ 2 

i - \u eh \ 2 



(37) 



A' reduces to the standard induced mass A = 2Ey/2GpN e in the limit where the heavy 
singlet neutrino decouples (\Uss\ — ► 1). The parameter 



B = 4EV2(G e N N e + G U N N U + G d N N d 



(3? 



describes the FCNC contributions from neutrino scattering off electrons and quarks in the 
sun or supernovae. It is convenient for later purposes to rewrite this as 



B = AEV2G F N e e(N n /N e ), 



(39) 



where 



e(x) = ^- \G% + 2G U N + G d N + (G U N + 2G d N )x 



(40) 



We have used the fact that the quark densities can be expressed in terms of N n and iV e 
(which equals the proton density in neutral matter). To compute B we set ^| = J2f e ^ c 
and obtain 



B = -2E\f2G F N n e ~ ±2Ey/2G F N n \U* h U t 



(h 



(41) 
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Note that due to the cancelation of V§ c and V^ c the parameter B is proportional to N n . 
Interestingly, this implies that the strongest possible effect arises from variation of ratio 

x = N n /N e . 



3 Constraints 



We now discuss experimental constraints on mixing between the ordinary neutrinos and 
SU(2)l singlets. A direct bound on {U^U^l comes from the KARMEN |19] experiment 
searching for "neutrino-oscillations" in the v e — 



channel. The upper bound on the 



transition probability P{y e 



yields 



w: h u, h \ 2 ^p{u e 



v. 



< 3.1 x KT 3 (90% c.l.; 



(42) 



This restricts \U* h Unh\ to be less than 0.06, which allows for rather large FCNCs. 



Much tighter bounds on mixing with a singlet neutrino are obtained from constraints on 
lepton universality, CKM unitarity, and the measured Z invisible decay width [pH[] . Rather 
than fitting the weak couplings of all known fermions we essentially update the analysis of 
Ref. [2I|, which considers the simpler possibility that only the neutrinos mix with exotic 
fermions. (|£/ a /i| 2 is equivalent to s 2 v in the notation of |[21]| .) The ratios of the leptonic 
couplings g e) g^ g T to the W are given by 




1- \U, 



ah 



u. 



a,b = e, /x, t. 



(43) 



bh 



The best test of e — fi universality of the W couplings comes from comparison of the 
rates for 7r — > fiu^ and n —* eu e , where the two most precise experiments J22| can be 
combined [^3|] to yield g^/ g e = 1.0012 ± 0.0016. The best tests of fi — r universality come 
from ALEPH measurements of leptonic 



and hadronic |25[ r decays. Averaging the 
two determinations of g T /gfj, (they are fully consistent with each other) yields g T /g^ = 
0.9943 ± 0.0065. From the leptonic decays one also obtains |24J g T /g e = 0.9946 ± 0.0064. 
The unitarity constraint for the first row of the CKM matrix is 



1=1,2,3 



(44) 



The most recent experimental determination P5j is X)i=i23 Wui\ 2 — 0.9972 ± 0.0013, 2a 
away from the SM value of unity. Finally, the ratio of the Z invisible decay width to the 
SM prediction is given by [|2"T| 



- 1 - Z—tinv 
Z— -Anv 



1 - 



\u eh \ 2 



Kh\ 2 



2\U Th \ 



(45) 
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Combining the most recent measurement of r z->mv [27] , obtained under the assumption 
of universal charged lepton couplings to the Z, with the SM prediction (for m t = 175.6 ± 
5.5 GeV) yields Y z ^ mv jYfi inv = 0.995 ± 0.004. 

We construct a \ 2 function with the experimental measurements discussed above, and 
derive bounds on the mixing parameters using the MINUIT package. Allowing for singlet 
mixing with all known neutrinos we obtain the 90% c.l. bounds 



\U eh \ 2 < 0.0049, \U Th \ 2 < 0.016, 



(46) 



and lUf.hl 2 = -0.0028 ± 0.0013, which is 2a away from the standard model but in the 
"wrong direction" for singlet mixing. This is due to the deviation of the CKM unitarity 
sum from unity. However, we note that this does not rule out z/^-singlet mixing since a 
discrepancy in this direction could be accounted for by mixing of the u and d quarks with 
SU{2) L singlets. We therefore present a second set of more conservative bounds (90% c.l.) 
in which all known neutrinos are allowed to mix with a singlet but the unitarity constraint 
has been eliminated: 



U eh \ 2 < 0.012, l^l 2 < 0.0096, \U Th \ 2 < 0.016. 



(47) 



We have not included the possibility of mixing between the charged leptons and exotics, 
nor have we taken into account correlations between ALEPH's determinations of g T /9fi 



and g T /g e in [23]. Nevertheless we can conclude that our parameter e should be less than 
about 1%, with the largest values possible corresponding to conversion of v e to v T . 



4 Solar Neutrino Oscillations 



Using the results of Section 2~5" it is straightforward to obtain the survival probability 
p8| , m, [18| for a v e that was produced in the solar center to be detected as an electron 
neutrino 



P N [v e -> u e ) = 
where the effective mixing is given by 

cos29n - 



^+ Q - P c ) cos 2ft cos 2(9 v. 



(4? 



(Acos2# - A') 



'(A cos 26 - A') 2 + (A sin 26 + B) 2 
The crossing probability P c is well approximated by ^9 



(49) 



P c ©(^4prod A res ) X 



exp[7r 7jv F(0)/2 


-exp 


'tt 1n F(6)/2 sin 2 26 


1 — exp 


'ti 1n F{6)/2 sin 2 20 





(50) 
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Here 7^ denotes the adiabaticity parameter which is given by [18 



A sin 2 29 

IN = 7TB oat a at 1 a \im\ x |l + 2e(x res )cot2g| 2 , (51) 

2Ecos29 (dN e /dx)/N e \ res 

where x res is the ratio N n /N p at the resonance. Using these results we have calculated 
the suppression rate for the three types of solar neutrino experiments as a function of A 
and sin 2 29. The calculations have been performed along similar lines to those described in 



Ref. fll8"H . We present the individual 95% c.l. contours (dashed for Homestake 0, dotted 
for the combined gallium experiments || |], and solid for Kamiokande ||) together with 
the combined allowed regions (shaded) for various positive (Fig. [TJ) and negative (Fig. ^) 
values of e. 

Although Z-induced FCNCs modify the individual contours at very small sin 2 29, the 
combined "small-angle solution" is (almost) unchanged for e = ±0.02 and the large-angle 
solution is not affected at all. In fact, even for larger e there are no dramatic changes 
to the standard small-angle MSW-solution. In light of the experimental constraints on e 
presented in Section ||] we conclude that Z-induced FCNCs must have negligible impact 
on solar neutrino oscillations. The cancelation of the NC potentials from electron and 
proton scattering leads to effects proportional to the neutron density N n , which in our sun 
is at most half of the charged particle densities N e = N p . This is unlike scenarios with 



flavor- changing neutrino interactions induced by new heavy particle exchange [30, 31, 32 



53 , M ], e.g., supersymmetric models without i?-parity, where the different contributions 
from scattering off quarks and electrons can add up constructively. 



5 Supernova Neutrino Oscillations 



Supernova explosions are intense neutrino sources. Due to the huge supernova densities 
these neutrinos can be resonantly converted for a large range of the parameters A and 
sin 2 29. The impact of neutrino oscillations has been discussed extensively in [^, |36, 
57 , [38]. Also it has been noted that neutrino oscillations might help to solve the shock 
reheating problem [[3£j 



In Refs. 40, 41 



resonant neutrino conversions were studied in 
the presence of supersymmetric i?-parity violating interactions and Z-induced FCNCs, 
but for massless neutrinos. In this Section we investigate how FCNCs due to mixing of 
the known neutrinos with singlets alter the MSW resonant conversion of massive neutrinos 
emerging from the neutrino-sphere of a supernova. The object we investigate is the survival 
probability P[y e ^e)- It has been shown in Ref. |36| that P[y e — > v e ) and P{p e —>■ u e ) are 
the only quantities that must be specified to determine how neutrino oscillations mix the 
fluxes. We assume that the neutrino parameters are such that matter-enhanced neutrino 
oscillations can occur for the neutrinos but not the anti-neutrinos. Moreover, we neglect 
neutrino-neutrino scattering which is justified for neutrino propagation outside the 
neutrino- sphere. For the point we wish to make it is sufficient to discuss the impact of 
FCNCs on the survival probability itself, so we do not convolute P{y e — > v e ) with the 
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predicted neutrino fluxes and cross sections in order to compute the experimental rates. 



(Note that from the neutrino data of SN1987a one cannot obtain a reasonable energy 
spectrum. Moreover, it might well be that NP processes, like Z-induced FCNCs, have 
an effect on the flux of produced neutrinos, which has not been considered in current 
supernova simulations [j43f .) 



As shown in Ref. JRJ], the following modifications of the solar MSW formalism are 
required in the case of supernova neutrinos: 

(a) The electron-density at the neutrino-sphere (N e ) pro< i ~ 10 35 cm -3 is larger than the solar 
core density (iV e ) 

core ~ 10 cm by about ten orders of magnitude. Since the supernova 
neutrinos are not much more energetic than solar neutrinos, an adiabatic threshold energy 



A cos 29 



2\^2GF(N e ) prod 



(52) 



of the order of a few MeV will scale the adiabatic band by a factor ~ 10 10 with respect 
to the solar case, to A max ~ 10 6 eV 2 . Above this value neutrinos cannot be resonantly 
converted, since they are produced below the resonance. Note that Ea is not changed by 
flavor- changing neutrino interactions 
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(b) Unlike the solar density profile, which decays (roughly) exponentially, the supernova 
density p(r) is predicted to decrease like 1/r 3 (r being the distance to the core outside 
the neutrino-sphere), but for the sake of generality we will just assume a power law (with 
a > 0): 



P(r) = P{R) 



R 



(53) 



For a density described by 
factor 7tv [defined in Eq. (| 
profile), but: 



3J) the scaling factor N' e /N e that appears in the adiabaticity 
1|)1 is not constant (like for an exponentially decaying density 



K 



-a: 



-a: 



R 



N e (r) 
N e (R) 



l/a 



(54) 



where we assume that the electron number density N e = pY e /m^ (m^ is the nucleon mass) 
is proportional to the mass density p (in fact Y e ~ 0.4 is constant to a good approximation 
outside the neutrino-sphere). To obtain 7^ we have to take N e in (^4j) at the resonance 



AT 



A cos 29 
2V2G F E' 



(55) 



resulting in 
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In 



sin 2 2^| 1 + 2t(x res )cot26\ 2 R f 2y/2G F EN e (R)' 



l/a 



2EA cos 29 

= I'— A -7- (cos 26) 
R r. 



a \ Acos26> / 
sin26> + 2e(x res )cos 29 | 2 



x- 



2a 



2 v / 2G F AT e ( J R) 



l/a 



(56) 



Thus, the effective non-adiabatic threshold energy is 



E 



NA 



A (cos 20)^ |sin2^ + 2e(x res )cos2^|^ x \2y/2G F N e (R) 



3.7 x 10 9 MeV 



/ A 



VeV z 



cos~ 2 29 I sin 29 + 2e(x res ) cos 201 



x 



4a 



(57) 



where for the last line we have taken typical supernova values, i.e., R = 10 7 cm for the 
radius of the neutrinosphere, p(R) = 10 12 g/cm 3 for the density at R and a = 3. From 
it is clear that the non-adiabatic band starts off at similar values of A m ,-, 



the sun, but has a different slope of —3/2 (for e(x r 
plane. 



^ 10 9 eV 2 as in 
,) = 0) in the (logarithmic) A — sin 2 29 



(c) Since the central supernova density is so huge, the "higher" e — r resonance almost 
always precedes the "lower" e — ft resonance. Thus the authors of Ref. have pointed 
out that a proper treatment of supernova neutrino oscillations should be done within a 
3-flavor formalism. Moreover they have noted in [33 that for small mixing angles the 
survival probability factorizes into 



PAv e -> u e ) x P h (u e 



(5? 



where P^h are the standard two-level survival probabilities for the lower and higher reso- 
nances, respectively. Since flavor- changing neutrino interactions become important when 
the mixing is small (i.e. tan 29 < |e|), eq. ([58]) is sufficient for our analysis. We will assume 
that A; and sin 2 29\ of the lower resonance are fixed by the standard MSW-solution to 
the solar neutrino problem (rather than fixing the ratios of A^/A^ and sin 2 29\j sin 2 26 \ at 
some arbitrary value) in order to obtain a two-dimensional plot. Then, at fixed energy E u , 
the survival probability P{y e — > u e ) is just a constant Pi multiplying Ph(E U) sin 26 h, A^, e) 
which we show in Fig. [3[ The plot exhibits the features of the supernova "MSW-triangle" 
that we discussed in (a) and (b): The adiabatic band appears at very large A max ~ 10 6 
eV 2 and the non-adiabatic band starts off at A min ~ 10~ 9 eV 2 for sin 2 26 = 1 extending to 
very small mixing sin 2 26 ~ 10~ 10 where it connects to the adiabatic band. 

The important consequence of this is that the supernova triangle is sensitive to even tiny 
FCNCs. One can see this clearly in Fig. f| and Fig. |5| where we show the survival probability 
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for various positive and negative e. The effect can be easily understood in terms of the 
adiabaticity parameter 7^. Without FCNCs 7 oc A^ 2 / 3 ) sin 2 26 for small vacuum mixing. 
Thus for smaller A (and fixed sin 2 26) the propagation is less adiabatic and there will be 
a minimal value A m j„ for each value of sin 2 26 where most of the neutrinos "cross-over" , 
resulting in a large survival probability. The non-adiabatic band is the contour defined 
by A m j„(sin 2 26) which separates the regions where the adiabatic conversion is efficient 
(above) and where it is not (below). As we have mentioned, this band is a straight line 
(with slope —3/2) if there are no FCNCs . However, if e 7^ then the adiabaticity parameter 
behaves like 7^ oc A*- 2 / 3 ^ sin2# + 2e| 2 for small vacuum mixing. Thus for sin 2$ <C |e|, A m j n 
is determined by e, the strength of the FCNCs , and not by the vacuum mixing sin 26 as can 
be seen in Fig. |] and Fig. [5]. Simply, in this regime the off-diagonal elements of the effective 
Hamiltonian in (|36D are dominated by the FCNC term rather than the mixing term. Note 
that for positive e (corresponding to negative e) the two competing contributions from 
FCNCs and mixing can cancel each other resulting in a vanishing 77V which implies a large 
survival probability around sin 26^ = — 2e (for a more detailed discussion see Ref. fl"8|l). 
The important result is that neutrino FCNC effects can be very significant for supernova 
neutrinos. 



6 Conclusions 

We have argued that neutrino singlets are almost unavoidable in any framework that at- 
tempts to solve neutrino anomalies by neutrino oscillations. The additional heavy neutrinos 
give rise to Z-induced FCNC interactions in the effective matter propagation matrix, and 
the question of whether or not they can be neglected is rather a quantitative than qual- 
itative one. We have worked out the resulting modifications to the MSW mechanism in 
order to study FCNC effects on matter-enhanced neutrino oscillations in the sun and in 
supernovae. We have found that while phenomenological constraints rule out significant 
changes in the solar neutrino MSW-solutions, the impact of FCNCs on the survival prob- 
ability of supernova neutrinos can be large, even for very small singlet components in the 
standard neutrino mass eigenstates. 

We conclude by asking whether values of e which are large enough to be of relevance to 
supernova neutrinos naturally occur in scenarios employing a neutrino see-saw mechanism. 
For example, this is easily seen not to be the case for the original see-saw matrix in eq. (§). 
The singlet-doublet mixing is of order mo/wiR so that e ~ m 2 D /m 2 R . This should be ^ 10~ 5 
for light-heavy mixing effects to be relevant for supernova neutrinos. But in this case the 
light neutrino masses would only be suppressed by a factor m^/m^ > 10~ 3 with respect 
to the Dirac masses. If the latter are reasonably large, e.g., of the order of the charged 
lepton masses, this would be inconsistent with the ultra-light neutrinos usually required to 
explain the solar neutrino problem (although strictly speaking only A has to be tiny). 

However, as already noted in the introduction, there are variants of the above scenario 
which can lead to much larger singlet admixtures. In a popular alternative, the Majorana 
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"double see-saw" [f45H , two singlets are added (per generation) leading to a 9 x 9 mass 
matrix of the form 



-A I' = j M T D AI n j . (59) 




For morris -C m R (mj = (detAfx) 1 ^ 3 ) singlet-doublet mixing is still of order mp/m^ 
but the light neutrino masses are of order ms^o/mii) 2 - Hence in this framework one can 
obtain ultra-light neutrinos while having a ratio mo/ fan > 10~ 3 that induces FCNCs that 
are significant for the MSW-effect in supernovae. Since the Dirac masses of the neutrinos 
have to be smaller than the weak scale, m^, it follows that as long as the mixing is 
~ m D /m R the Z-induced FCNCs are potentially relevant to our discussion if m R < 100 
TeV. 

Finally, vectorlike pairs of SU{2) L singlet quarks and leptons with large masses m v are 
often introduced in order to suppress known quark and charged lepton masses relative to the 



weak scale via a generalized "Dirac see-saw" [46], leading to left-handed singlet components 



in the ordinary charged fermion mass eigenstates of order m wea k/mv- If vectorlike pairs of 
neutral singlets are also present in such a scenario the known neutrinos could be expected 
to mix with the left-handed singlets at same order as the charged fermions, in addition 
to mixing with right-handed singlets responsible for an ultra-light Majorana mass see- 
saw. (In the absence of a Majorana see-saw the vectorlike singlets would typically lead 
to neutrino Dirac masses which are of same order as the charged lepton masses.) As 
in the above example, singlet mixing at the level of interest for supernova oscillations, 
i.e., m weak /m v > 10~ 3 , would correspond to a New Physics mass scale m v < 100 TeV. 



The right-handed neutrino Majorana mass scale m R could be arbitrarily large f47fl , thus 
allowing for ultra-light neutrino masses which are consistent with the MSW solution for 
solar neutrinos. 



Note added: When this work was near completion we learned of another paper |48[ that 
analyzed the effects of FCNCs on supernova oscillations in the context of supersymmetric 
models with broken .R-parity, arriving at effects of similar magnitude to those presented in 
our analysis. 
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Figure 1: The combined allowed regions for the solar neutrino experiments with e > 
as indicated above each plot. The dotted contours correspond to the combined gallium 
experiments, the dashed contours to the Homestake experiment and the solid contours to 
the Kamiokande experiment. The shaded areas indicate the 95% c.l. combined allowed 
regions in the sin 2 29 — A plane. 

19 



A [eV 2 ] 



e=0 




w w 10 10 10 




10 



£=-0.02 




£=-0.05 




10 10 




sin 2 29 

Figure 2: The combined allowed regions for the solar neutrino experiments with e < 
as indicated above each plot. The dotted contours correspond to the combined gallium 
experiments, the dashed contours to the Homestake experiment and the solid contours to 
the Kamiokande experiment. The shaded areas indicate the 95% c.l. combined allowed 
regions in the sin 2 29 — A plane. 
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Figure 3: The MSW-contours for supernova neutrinos with e = at one discrete energy 
(E u = 10 MeV). The shading indicates the value of the survival probability PyXye v e) 
in the sin 2 29 — A plane: White corresponds to 0.9 < Ph < 1.0 and the darkest area 
corresponds to 0.0 < Ph < 0.1. 
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Figure 4: The MSW-contours for supernova neutrinos with e > (indicated above each 
plot) at one discrete energy [E v = 10 MeV). The shading indicates the value of the survival 
probability Ph(y e — > v e ) in the sin 2 29 — A plane: White corresponds to 0.9 < Ph < 1.0 
and the darkest area corresponds to 0.0 < Ph < 0.1. 
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Figure 5: The MSW-contours for supernova neutrinos with e < (indicated above each 
plot) at one discrete energy [E v = 10 MeV). The shading indicates the value of the survival 
probability Ph(y e — > v e ) in the sin 2 29 — A plane: White corresponds to 0.9 < Ph < 1.0 
and the darkest area corresponds to 0.0 < Ph < 0.1. 
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